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Abstract. In this paper, we consider the problem of choosing disks (that 
we can think of as corresponding to wireless sensors) so that given a set 
of input points in the plane, there exists no path between any pair of 
these points that is not intercepted by some disk. We try to achieve this 
separation using a minimum number of a given set of unit disks. We show 
that a constant factor approximation to this problem can be found in 
polynomial time using a greedy algorithm. To the best of our knowledge 
we are the first to study this optimization problem. 



1 Introduction 

Wireless sensors are being extensively used in applications to provide barriers 
as a defense mechanism against intruders at important buildings, estates, na- 
tional borders etc. Monitoring the area of interest by this type of coverage is 
called barrier coverage [11]. Such sensors are also being used to detect and track 
moving objects such as animals in national parks, enemies in a battlefield, forest 
fires, crop diseases etc. In such applications it might be prohibitively expensive 
to attain blanket coverage but sufficient to ensure that the object under consid- 
eration cannot travel too far before it is detected. Such a coverage is called trap 
coverage [3, 14]. 

Inspired by such applications, we consider the problem of isolating a set of 
points by a minimum-size subset of a given subset of unit radius disks. A unit 
disk crudely models the region sensed by a sensor, and the work reported here 
readily generalizes to disks of arbitrary, different radii. 

Problem Formulation. The input to our problem is a set I of n unit disks, and 
a set P of k points such that I separates P, that is, for any two points p,q e P, 
every path between p and q intersects at least one disk in I. The goal is to 
find a minimum cardinality subset of I that separates P. See Figure 1 for an 
illustration of this notion of separation. 

There has been a lot of recent interest on geometric variants of well-known 
NP-hard combinatorial optimization problems, and our work should be seen in 
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this context. For several variants of the geometric set cover problem, for example, 
approximation algorithms have been designed [7, 2, 12] that improve upon the 
best guarantees for the combinatorial set cover problem. For the problem of 
covering points by the smallest subset of a given set of unit disks, we have 
approximation algorithms that guarantee an O(l) approximation and even a 
PTAS [4, 12]. These results hold even for disks of arbitrary radii. Our problem 
can be viewed as a set cover problem where the elements that need to be covered 
are not points, but paths. However, known results only imply a trivial O(n) 
approximation when viewed through this set cover lens. 

Another example of a problem that has received such attention is the inde- 
pendent set problem. For many geometric variants [5, 6, 9], approximation ratios 
that are better than that for the combinatorial case are known. 

Our problem is similar to the node multi-terminal cut problem in graphs 
[10]. Here, we are given a graph G = (V,E) with costs on the vertices and a 
subset U C V of k vertices, and our goal is to compute a minimum cost subset 
of vertices whose removal disconnects every pair of vertices in U . This problem 
admits a poly-time algorithm that guarantees an O(l) approximation. We note 
however that the problem we consider does not seem to be a special case of the 
multi-terminal cut problem. 

Contribution and Related Work. Our main result is a polynomial time algorithm 
that guarantees an O(l) approximation for the problem. To the best of our 
knowledge, this is the first non-trivial approximation algorithm for this problem. 
Our algorithm is simple and combinatorial and is in fact a greedy algorithm. 
We first present an O(l) approximation for the following two-point separation 
problem. We are given a set of unit disks G, and two points s and t, and we wish 
to find the smallest subset BCGso that B separates s and t. 

Our greedy algorithm to the overall problem applies the two-point separation 
algorithm to find the cheapest subset B of / that separates some pair of points 
in P. Suppose that P is partitioned into sets Pi,p2, ■ . ■ ,P T where each P, is 
the subset of points in the same "face" with respect to B. The algorithm then 
recursively finds a separator for each of the Pi, and returns the union of these 
and B. 
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The analysis to show that this algorithm has the 0(1) approximation guar- 
antee relies on the combinatorial complexity of the boundary of the union of 
disks. It uses a subtle and global argument to bound the total size of all the 
separators B computed in each of the recursive calls. 3 

Our approximation algorithm for the two-point separation problem, which 
is a subroutine we use in the overall algorithm, is similar to fast algorithms 
for finding minimum s-t cuts in undirected planar graphs, see for example [13]. 
Our overall greedy algorithm has some resemblance to the algorithm of Erickson 
and Har-Peled [8] employed in the context of approximating the minimum cut 
graph of a polyhedral manifold. The details of the our algorithm and the analysis, 
however, are quite different from these papers since we do not have an embedded 
graph but rather a system of unit disks. Sankararaman et al. [14] investigate a 
notion of coverage which they call weak coverage. Given a region 1Z of interest 
(which they take to be a square in the plane) and a set / of unit disks (sensors) , 
the region is said to be fc-weakly covered if each connected component of 1Z — 
U deI d has diameter at most k. They consider the situation when a given set / 
of unit disks completely covers 1Z, and address the problem of partitioning / into 
as many subsets as possible so that 1Z is k- weakly covered by every subset. Their 
work differs in flavor from ours mainly due to the assumption that / completely 
covers 1Z. 

Organization. In Section 2, we discuss standard notions we require, and then 
reduce our problem to the case where none of the points in input P are contained 
in any of the input disks. In Section 3, we present our approximation algorithm 
for separating two points. In Section 4, we describe our main result, the constant 
factor approximation algorithm for separating P. We conclude in Section 5 with 
some remarks. 

2 Preliminaries 

We will refer to the standard notions of vertices, edges, and faces in arrangements 
of circles [1]. In particular, for a set R of m disks, we are interested in the faces 
in the complement of the union of the disks in R. These are the connected 
components of the set 3? 2 — \J d&R d. We also need the combinatorial result that 
the number of these faces is 0(m). Furthermore, the total number of vertices and 
edges on the boundaries of all these faces, that is, the combinatorial complexity 
of the boundary of the union of disks in R, is 0(m) [1]. We make standard 
general position assumptions about the input set / of disks in this article. This 
helps simplify the exposition and is without loss of generality. 

3 In the earlier version of this paper, a similar algorithm was analyzed in a more 
"local" fashion. The basic observation was that the very first separator B that is 
computed has size 0(\OPT\/k), where OPT is the optimal solution for the problem. 
Subsequent separators computed in the recursive calls may be more expensive, but it 
was shown that the overall size is 0((log k) ■ \OPT\). In contrast, the present analysis 
does not try to bound the size of the individual separators, but just the sum of their 
sizes. As a consequence, the analysis also turns out to be technically simpler. 
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Lemma 1. Let R be a set of disks in the plane, and Q a set of points so that 
(a) no point from Q is contained in any disk from R, and (b) no face in the 
complement of the union of the disks in R contains more than one point of Q. 
Then \R\ = Q{\Q\). 

Proof. The number of faces in the in the complement of the union of the disks 
in R is 0(\R\). □ 

Covering vs. Separating. The input to our problem is a set / of n unit disks, 
and P a set of k points such that / separates P. Let P c C P denote those points 
contained in some disk in /; and P s denote the remaining points. We compute 
an a-approximation to the smallest subset of / that covers P c using a traditional 
set-cover algorithm; there are several poly-time algorithms that guarantee that 
a = 0(1). We compute a f3- approximation to the smallest subset of / that 
separates P s , using the algorithm developed in the rest of this article. We argue 
below that the combination of the two solutions is an 0(a + /?) approximation 
to the smallest subset of I that separates P. 

Let OPT C / denote an optimal subset that separates P. Suppose that OPT 
covers k\ of the points in P c and let k 2 = \Pc\ — k\. By Lemma 1, |OPT| = Q(k2). 

Now, by picking one disk to cover each of the k 2 points of P c not covered by 
OPT, we see that there is a cover of P c of size at most \OPT\ + k 2 = 0(\OPT\). 
Thus, our a-approximation has size 0(a) ■ \OPT\. Since OPT also separates P s , 
our /3-approximation has size 0((3) ■ \ OPT\. Thus the combined solution has size 
0(a + /3) ■ \OPT\. 

In the rest of the article, we abuse notation and assume that no point in the 
input set P is contained in any disk in /, and describe a poly-time algorithm 
that computes an 0(l)-approximation to the optimal subset of / that separates 
P. 



3 Separating Two Points 

Let s and t be two points in the plane, and G a set of disks such that no disk in 
G contains either s or t, but G separates s and t. See Figure 2. Our goal is to 
find the smallest cardinality subset B of G that separates s and t. We describe 
below a polynomial time algorithm that returns a constant factor approximation 
to this problem. 

Without loss of generality, we may assume that the intersection graph of G is 
connected. (Otherwise, we apply the algorithm to each connected component for 
which the disks in the component separate s and t. We return the best solution 
obtained.) Let f s and ft denote the faces containing s and t, respectively, in 
the arrangement of G. We augment the intersection graph of G with vertices 
corresponding to s and t, and add an edge from s to each disk that contributes 
an edge to the boundary of the face f s , and an edge from t to each disk that 
contributes an edge to the boundary of the face f t . We assign a cost of to s, t, 
and a cost of 1 to each disk in G. We then find the shortest path from s to t in 
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this graph, where the length of a path is the number of the vertices on it that 
correspond to disks in G. Let a denote the sequence of disks on this shortest 
path. Note that any two disks that are not consecutive in a do not intersect. 

Using a, we compute a path ir in the plane, as described below, from s to t 
so that (a) there are points s' and t' on 7r so that the portion of 7r from s to s' 
is in / s , and the portion from t' to t is in f t ; (b) every point on it from s' to t' 
is contained in some disk from a; (c) the intersection of tt with each disk in a is 
connected. See Figure 2. 

Suppose that the sequence of disks in a is d\, . . . , d\ a \ . Let s' (resp. t') be a 
point in d\ (resp. d a ) that lies on the boundary of f 8 (resp. f t ). For 1 < i < \a\ — 1, 
choose to denote an arbitrary point in the intersection of di and d^i . The path 
7r is constructed as follows: Take an arbitrary path from s to s' that lies within 
/ s , followed by the line segments s'x\,x±X2, • ■ • , x\cr\-2 x \<r\-i> x \a\ -it' i followed 
by an arbitrary path from t' to t that lies within ft . 

Properties (a) and (b) hold for n by construction. Property (c) is seen to 
follow from the fact that disks that are not consecutive in a do not overlap. 

Notice that 7r "cuts" each disk in a into two pieces. (Formally, the removal of 
7T from any disk in a yields two connected sets.) The path 7r may also intersect 
other disks and cut them into two or more pieces, and we refer to these pieces 
as disk pieces. For a disk that 7r does not intersect, there is only one disk piece, 
which is the disk itself. 

We consider the intersection graph H of the disk pieces that come from disks 
in G. Observe that a disk piece does not have points on ir, since tt is removed; 
so two disk pieces intersecting means there is a point outside n that lies in both 
of them. In this graph, each disk piece has a cost of 1. 

In this graph H, we compute, for each disk d G a, the shortest path between 
the two pieces corresponding to d. Suppose d! € o yields the overall shortest path 
a'; let D denote the set of disks that contribute a disk piece to this shortest path. 
Our algorithm returns D as its computed solution. See Figure 3. 

We note that D separates s and t - in particular, the union of the disk pieces 
in a' and the set 7T R d' contains a cycle in the plane that intersects the path w 
between s and t exactly once. 
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Fig. 3. This figure continues with the example of Figure 2. The disks with bold bound- 
ary are the set D computed by our algorithm. The only disk from G with bold boundary 
has two disk pieces, and the shortest path between them in graph H yields D. 



3.1 Bounding the Size of the Output 



Let B* denote the smallest subset of G that separates s and t. We will show 
that \D\ = 0(\B*\). Let /* denote the face containing s in the arrangement of 
B*. Due to the optimality of B* , we may assume that the boundary of /* has 
only one component. Let a (resp. b) denote the first (resp. last) point on path 
7r where tt leaves /*. It is possible that a = b. Wc find a minimum cardinality 
contiguous subsequence <f of a that contains the subpath of 7r from a to b; let 
d a and c4 denote the first and last disks in a. See Figure 4. 




(a) (b) 

Fig. 4. The shaded disks are in _B*.(a) This figure shows points a and b where n leaves 
/* for the first and last time, respectively, (b) This figure shows the face / containing 
s in the arrangement with B* U a 
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We claim that \a\ < \B* \ + 2; if this inequality does not hold, then we obtain 
a contradiction to the optimality of a by replacing the disks in the a\ {d a ,db} 
by P*. 

Consider the face / containing s in the arrangement with B* U a. Each 
edge that bounds this face comes from a single disk piece, except for one edge 
corresponding to d a that may come from two disk pieces. (This follows from the 
fact that the portion of ir between a and b is covered by the disks in a.) These 
disk pieces induce a path in H in between the two pieces from d a , and their 
cost therefore upper bounds the cost of D. We may bound the cost of these disk 
pieces by the number of edges on the boundary of / (with respect to B* U a). 
The number of such edges is 0(\B* \ + \a\) = 0(\B*\). 

Theorem 1. Let s and t be two points in the plane, and G a set of disks such 
that no disk in G contains either s or t, but G separates s and t. There is a 
polynomial time algorithm that takes such G, s, and t as input, and outputs a 
subset B C G that separates s and t; the size of B is at most a multiplicative 
constant of the size of the smallest subset B* C G that separates s and t. 

4 Separating Multiple Points 

We now present a polynomial time algorithm that yields an 0(1) approximation 
to the problem of finding a minimum subset of I that separates the set P of 
points. The algorithm is obtained by calling recSep(P), where recSep(Q), for 
an y Q Q P is the following recursive procedure: 

1. If \Q\ < 1, return 0. 

2. For every pair of points s,t 6 Q, invoke the algorithm of Theorem 1 (with 
G •<— I) to find a subset B s t C I such that B s t separates s and t. 

3. Let B denote the minimum size subset B Sit over all pairs s and t considered. 

4. Consider the partition of Q into subsets so that each subset corresponds to 
points in the same face (with respect to B). Suppose Qi, . . . ,Q T are the 
subsets in this partition. Note that r > 2, since B separates some pair of 
points in Q. 

5. Return B U Uj=i recSep(Qj). 

Clearly, recSep(P) yields a separator for P. To bound the size of this sepa- 
rator, let us define a set Q that contains as its element any Q C P such that 
\Q\ > 2 and recSep(Q) is called somewhere within the call to recSep(P). For any 
Q £ Q, define Bq to be the set B that is computed in the body of the call to 
recSep(Q). Notice that recSep(P) returns Uq 6 qPq. 

Now we "charge" each such Bq to an arbitrary point within pq £ Q in such 
a way that no point in P is charged more than once. A moment's thought reveals 
that this is indeed possible. (In a tree where each interval node has degree at 
least 2, the number of leaves is greater than the number of internal nodes.) 

Let OPT denote the optimal separator for P and let Fq C OPT denote the 
disks that contribute to the boundary of the face (in the arrangement of OPT) 
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containing pq. We claim that \Bq \ = 0(\Fq\); indeed Fq separates pq £ Q from 
any point in P, and thus any point in Q. Thus for any t £ Q \ {pq}, we have 
\Bq\ <B PQ<t = 0{\F Q \), 
We thus have 

(J \B Q \<J20(\F Q \) = 0(\OPT\), 
QeQ QeQ 

where the last equality follows from union complexity. 
We have derived the main result of this paper: 

Theorem 2. Let I be a set of n unit disks and P a set of k points such that I 
separates P. There is a polynomial time algorithm that takes as input such I and 
P, and returns a subset O C / of disks that also separates P , with the guarantee 
that \0\ is within a multiplicative 0(1) of the smallest subset of I that separates 
P. 

5 Conclusions 

Wc have a presented an 0(l)-approximation algorithm for finding the minimum 
subset of a given set / of disks that separates a given set of points P. One way 
to understand our contribution is as follows. Suppose we had at our disposal 
an efficient algorithm that optimally separates a single point p £ P from every 
other point in P. Then applying this algorithm for each point in P, we get a 
separator for P. That the size of this separator is within 0(1) of the optimal is 
an easy consequence of union complexity. However, we only have at our disposal 
an efficient algorithm for a weaker task: that of approximately separating two 
given points in P. What we have shown is that even this suffices for the task of 
obtaining an O(l) approximation to the overall problem. 

It is easy to see that our algorithm and the approximation guarantee gen- 
eralize, for example, to the case when the disks have arbitrary and different 
radii. 

Acknowledgements. We thank Alon Efrat for discussions that led to the for- 
mulation of the problem, and Sariel Har-Peled for discussions that led to the 
algorithm described here. 
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